Conversion of hole states by acoustic solitons 
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The hole states in the valence band of a large class of semiconductors are degenerate in the 
projections of angular momentum. Here we show that the switching of a hole between the states 
can efficiently be realized by acoustic solitons. The microscopic mechanism of such a state conversion 
is related to the valence band splitting by local elastic strain. The conversion is studied here for 
heavy holes localized at shallow and deep acceptors in silicon quantum wells. 

PACS numbers: 78.20.Hp, 63.20.K-, 78.67.De 



The strain induced effects have been proved to be an 
efficient tool to study basic properties of semiconduc- 
tor structures yielding information on their point-group 
symmetry, band parameters, spin-orbit coupling, etc.— 
During the last decade the possibility emerged to study 
such effects in nanostructures on the extremely short time 
scale of picosecond range. This could be done employ- 
ing acoustic solitary waves (solitons) i^i^i^ The acoustic 
soliton is a pulse of strong nonlinear elastic strain propa- 
gating along the crystal with the velocity larger than the 
sound speedi^ Due to rather high values of deformation- 
potential constants in semiconductors, such a local strain 
efficiently modifies the band structure and thereby inter- 
acts with charge carriers. The important point is that, 
unlike linear acoustic waves, the soliton-type strain does 
not change its sign within the pulsei^ii Therefore, after 
the soliton has passed through the sample, the electron 
state of a system may differ from the initial one. 

In this paper we show that, in semiconductor struc- 
tures where the electron or hole states are degenerate 
in the absence of strain, acoustic solitons can cause a 
transition of carriers between these quantum-mechanical 
states. We study such a state conversion for heavy holes 
localized at acceptors in quantum wells (QWs). The 
ground state of localized holes in most of the cubic semi- 
conductors is degenerate in the projection of the angular 
momentum because of the complex structure of the va- 
lence band. We show that the propagation of the acoustic 
soliton of a certain amplitude through the area of hole 
localization changes the projection of the hole angular 
momentum. The effect of state conversion, particularly 
between two easily distinguished states, is of interest not 
only from fundamental point of view but can also be uti- 
lized for the information processing and storage. 

The acoustic soliton represents a perturbation of 
purely mechanical origin, therefore, in the first approx- 
imation it does not interact with the carrier spin. Ac- 
cordingly, it is convenient to consider the effect of state 
conversion for semiconductor structures with negligible 
spin-orbit interaction (such as Si, SiC, etc). Below we 
focus on silicon-based quantum wells although the main 
conclusions can be generalized to other systems. We as- 



sume that in bulk material the valence-band states at the 
center of the Brillouin zone belong to the representations 
r'js or and denote the basis functions as X, Y, and 
Z. Due to quantum confinement, the ground heavy-hole 
states in a structure grown along the z axis are formed 
from the Bloch amplitudes X and Yi^ The corresponding 
wave functions have the form 



*(r) a{p)u{z)X + (i{p)u{z)Y 



(1) 



where a{p) and /3(p) arc smooth envelopes in the QW 
plane, which can be combined into a two-component col- 
umn il^ip) ~ [o^ip] 1 P{p)V : P = {x,y) is the in-plane 
coordinate, u(z) is the function of size quantization, x, 
y, z are the cubic axes, and the spin index is omitted. 

The envelope functions V'(p) of a hole localized at an 
acceptor in a narrow quantum well can be found by solv- 
ing the matrix Schrodinger equation Hip{p) = Eip{p), 
where H is the effective Hamiltonian assuming in the 
spherical approximation the form 
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Here k = (kx, ky) is the momentum operator divided by 
the reduced Planck constant h, 7 and 7' are the band 
parameters which are expressed in terms of the param- 
eters L and M (see Ref. [l|) via 7 = (L + M)/2 and 
7' = (L — M)/2, U{p) is the attractive potential of 
the acceptor, and / is the unit matrix 2x2. We as- 
sume the acceptor potential U{p) to be isotropic. Then, 
the hole states are described by the quantum number n 
{n = 0, 1, 2, . . .) and the projection of the orbital angular 
momentum onto the z axis j (j = 0, ±1, ±2, . . .). The 
corresponding wave functions ipnj{p) in the polar coor- 
dinates p = (/9, ifi) are given by 



i^njip) 



■Hnjip) + £,nj{p)' 

iVnjip) - i&Ljip)e^"^ 



(3) 



where r]nj (p) and (p) are solutions of the equation set 
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Vnjip)^[U{p)-Enj]^nj{p), (4) 
Vnj{p) = [U{p) - Enj]'nnj{p) , 



Enj is the energy of the state ipnj (p) measured from the 
subband edge. We note that the states {n,j) and (n, —j) 
have the same energy and their envelope functions are 
related by ■ij)„-j{p) = ■ip^.j{p). From Eq. ^ it follows 
that only the states with j = ±1 contain envelope func- 
tions of s-typc which are nonzero at the acceptor. There- 
fore, the ground state of localized holes is formed with 
the functions ^o.±i(p)- Below we denote them as ijj±{p), 
respectively. 

In the absence of strain, the ground state is degener- 
ate in the projection of angular momentum and the hole 
wave function represents a superposition of V'+(p) ^^id 
■0- (p) . The acoustic soliton lifts the degeneracy causing 
the transition of a hole between the above states. 

The strain effect on the heavy-hole subband is de- 
scribed by the effective Hamiltoniar^i 
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lUxx + T^T-U 
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(5) 



where m, and n are the deformation potential con- 
stants, Uq/3 are the strain tensor components. We con- 
sider that the bulk acoustic soliton propagates in the x 
direction along the quantum well plane inducing time- 
dependent component Uxx{x,t) of the strain tensor, see 
inset to Fig. [1] The strain pulse is assumed weak enough 
not to cause ionization of the localized hole or its tran- 
sition to excited states. Then, according to the pertur- 
bation theory, the hole wave function ip{p^t) can be ex- 
panded over the non-disturbed states V'±(p) as follows 



V'(p,<) - c+(t)v+(p) + c_(t)v_(p). 



(6) 



where the coefficients Cj{t) {j = ±) satisfy the coupled 
equations 

ih^^ = V++{t)c+{t) + V+.{t)c^{t), (7) 



dt 

dt 



= l/_+(i) c+(i) + ^__(i)c_(i) 



and Vjj' (t) ~ J 'ip'j{p)Vipj' {p)dp are the matrix elements 
of the perturbation ([5]). The specific form of the func- 
tions ([3]) and the perturbation V oc Uxxix,t) leads to the 
relations V++{t) = V--{t) and V+-{t) = V-+{t), which 
makes Eqs. (O easily solvable. 

The solution of Eqs. ([7]) assumes the form 

c+{t) = [a+cos$(t)-m_sin$(t)]c-*®(*) , (8) 
c-{t) = [a_ cos$(t) - m+sin$(i)] e"*®^*) , 



where a+ = c+(— oo) and a_ = c_(— cx)) are coefficients 
of the initial state at t — — oo, (f>(t) and 0(t) are phases 
given by 



v+^{t')dt' , e{t) 



V++{t')dt' . (9) 



It follows from Eqs. ([8]) that it is the phase ^{t) that 
describes the hole state evolution while Q{t) constitutes 
the common phase factor. In the final state, i.e., at t = 
+00 when the soliton has completely passed through the 
area of hole localization, the phase shift has the form 
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(10) 



Since the average value of the strain Uxx is non-zero 
within the soliton, <!>/ ^ 0, and the final hole state differs 
from the initial one. The most efficient conversion of the 
hole state occurs if the soliton amplitude is so that 



|$/| =7r(?i + l/2), 71 = 0,1,2,. 



(11) 



In this particular case, the soliton-hole interaction results 
in the complete switching of the hole angular momentum: 
The initial state ip+{p) is converted into the state ip-{p) 
and vice versa, see Eq. ([8]). 

To elaborate the effect of state switching we con- 
sider two models of the acceptor: (i) Coulomb potential 
U{p) = — e2/(ep), where e is the elementary charge and 
e is the dielectric constant, and (ii) zero-radius potcn- 
tialiiSiiii In the first case, the ground hole states and the 
localization energy ii^o > are calculated numerically by 
solving Eqs. The latter approach allows us to find 
the states analytically for a given Eq assuming that the 
range of attractive potential is much shorter than the 
hole localization length. Within this approach, solution 
of Eq. ([D) for the ground state has the form 
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where Kq{x) and if 2 (a;) are the Macdonald fmictions, the 
radii pi and p2 are given by pi = a/-(7 + j')/Eo, p2 = 
— ■y')/Eo, and it is assumed that 7 + 7', 7 — 7' < 0. 
Acoustic solitons in crystals arc typically approxi- 
mated by solutions of the Korteweg-de Vrics (KdV) wave 
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FIG. 1: (Color online) Dependence of the phase shift $/ 
on the strain amplitude uq calculated numerically for the 
Coulomb potential (squares) and plotted after Eq. (|14p for the 
zero-radius potential (solid curve). Dashed line corresponds 
to $/ = 7r/2 which is optimal for the state conversion. Inset 
sketches the propagation of strain soliton along the quantum 
well. 



equation^ or the doubly dispersive equation (DDE)i^ 
Since the particular soliton shape is not crucial for our 
results, wc consider the KdV soliton so that u^x has the 
form 



5(2;, t) = uq cosh 



vt 



L 



(13) 



where uq is the strain amplitude, v and L = d/y/u^ are 
the soliton velocity and size, respectively, and d is a ma- 
terial constant. In this case j^'^ u^xix, i)dt = Id^fu^jv. 

Figure [T] shows dependence of the phase shift <t>/ on the 
strain amplitude uq. Squares correspond to the numeri- 
cal calculation of <t>/ given by Eq. (|10p for the Coulomb 
localizing potential, solid curve represents the analytical 
dependence 
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derived in the zero-radius-potcntial approach by inte- 
grating the wave function (|12|) . The band parame- 
ters of silicon used in the calculation arc as follows: 
7 = -4.65 h/{2mo), 7' = -1.13 h/{2mo)^ where toq is 
free electron mass, the deformation-potential constants 
I = —4.9 eV and m = —1.5 eVr^ The soliton parameters 
d = 1.5 A, w = 0.9-10^ cm/s, and the strain amplitude uq 
in the range 10~^-^10~'^ are chosen, which corresponds to 
experimental data on acoustic phonon pulses in silicon^i^ 
From Fig. [T] it follows that the phase shift = 7r/2 
optimal for the state conversion is achieved in Si-based 
structures at a moderate strain amplitude uq w 1.8-10"'^. 
It proves that the hole states can be efficiently manipu- 
lated by acoustic solitons. Moreover, both Coulomb and 
zero-radius potentials lead to the same quantitative re- 
sult indicating that $ / weakly depends on the form of lo- 
calizing potential for silicon band parameters. This can 
be attributed to the fact that the envelope function of 
the ground heavy-hole state is mainly of s-type even at 
7'/7 « 0.24. Therefore, the integral / rj^ i{p)dp deter- 



mining the conversion efficiency [see Eq. (jlOp ] is approx- 
imately equal to 1 in accordance with the wave function 
normalization and independent of the explicit form of 

For the strain amplitude uq = 1.8 • 10^''^ optimal for 
the state conversion, the soliton length L is estimated 
as 110 A that is approximately 7 times larger than the 
radius of hole localization aq = e^/ (ei?o) ~ 16 A in silicon 
quantum wells. In the approximation of L 3> ao, the time 
evolution of the phase <I>(t) for the soliton shape has 
the form 



$(i) = $ 



tanh (vt/L) + 1 



f- 



(15) 



The soliton-induced evolution of a hole state can fruit- 
fully be considered as a precession of pseudospin or a 
trajectory on the Bloch sphere. In this approach^i^ any 
superposition ([6]) of two states, ifj+ip) and ^p-{p), is at- 
tributed to a unit vector S or a point {Sx, Sy, Sz) on the 
sphere, see Fig. [21 The components of S are given by 



(16) 



where cr is the vector of Pauli matrices and x is the 
spinor composed of the coefficients c+(i) and c_(i), 
X = [c+{t),c^{t)]'^ . In particular, the pure states 
^+{p) and ip-{p) correspond to the polar points (0, 0, 1) 
and (0,0,-1), respectively, while the states ['0+(p) + 
e^°'ifi^{p)]/\/2 with arbitrary phase a are projected onto 
the Bloch sphere equator. 

During the soliton-hole interaction the coefficients c± 
vary in time, see Eq. ([5]), which corresponds to motion 
of the point {Sx, Sy, Sz) on the Bloch sphere. Equation 
describing the time evolution of S can be derived from 
Eq. (O. It gives 



- = m) X s] 



(17) 



where n{t) ~ (l/^)X]j/ ^jj'i^) ^j' j- ^'^'^ ^ given initial 
state Sq and 0,{t), Eq. pT]) allows one to calculate the 
trajectory S{t) and restore the wave function. Shown in 
Fig. [2] are examples of such trajectories under a soliton- 
induced uniaxial strain. The strain Uxxix,t) produced 
by soliton propagating along the x axis corresponds to 
n{t) II X and, therefore, causes S to rotate around the 
X axis (Fig. [2^) while the strain Uyy{y,t) corresponds to 
r2(t) II y leading to the rotation of S around the y axis 
(Fig. Wp)- In the case of complete conversion of the hole 
state from ip-\-{p) into ijj^{p), the trajectories pass from 
one pole of the Bloch sphere to the other. 

The time of switching the hole states is given by 
T = 2(L -|- ao)/v. It is of nanosecond scale and much less 
than the lifetime of localized holes at low temperatures as 
well as the soliton lifetime. It suggests that the effect can 
be studied experimentally by means of time- and space- 
resolved optical spectroscopy. Indeed, the states ip+{p) 
and ip^{p) are characterized by projections -|-1 and —1 of 
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FIG. 2: (Color online) Trajectories on the Bloch sphere cor- 
responding to the time evolution of the hole state caused by 
acoustic soliton propagating along (a) x and (b) y axes. 



the orbital angular momentum, therefore, optical pump- 
ing by circularly polarized light leads to a predominant 



population of one of the states. The hole angular mo- 
mentum can be registered in turn by analyzing the po- 
larization of recombination radiation or the polarization 
change of a probe pulse (see, e.g., Ref. The soliton 

propagation through the area of hole localization reverses 
the projection of hole angular momentum leading to a 
change in optical response. 

To summarize, we have demonstrated the possibility 
of switching the quantum-mechanical state of a localized 
hole by an acoustic soliton. The amplitude of the strain 
soliton required for switching the hole states in silicon- 
based quantum wells is found to be « 2 • 10~^ which 
corresponds to strain pulses studied experimentally. 
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